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5!!^ . Abstract 

^ ■ 
o ' We propose a model for scattering in a flat resonator with a thin 

antenna. The results are applied to rectangular microwave cavities. 

G I We compute the resonance spacing distribution and show that it agrees 

^ ' well with experimental data provided the antenna radius is much 

O^' smaller than wavelengths of the resonance wavefunctions. 

> 

X 

>y , The aim of the present letter is twofold. First of all, we want to present a 

method to treat scattering in systems consisting of components of different 
dimensionality. To illustrate how it works, we are going then to apply it to 
the case of a flat microwave cavity to which a thin antenna is attached. We 
compute the resonance spacing distribution in this system and show that it 
reproduces experimental data — within certain range of energies but without 
introducing any free parameters. 

Many recent studies involve an analysis, both theoretical and experi- 
mental, of spectral and transport properties of systems with a complicated 
geometry: let us recall various microwave resonators [0, H, §, ^, ^ or conduc- 
tance fluctuations in quantum dots — see [^, 0, H, |^ and references therein. 
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Although physically different such systems share many common properties, 
in particular, chaotic spectral behavior for certain "resonator" shapes. For 
the sake of definiteness, we shall speak here only about the electromagnetic 
case considering the experimental setup of the work mentioned above, i.e., a 
fiat microwave cavity which may be regarded as two-dimensional coupled to 
one or more antennas that supply power compensating radiative loss. 

Already the very first experiments 0] revealed a discrepancy with theo- 
retical predictions: the expected Poissonian eigenvalue-spacing distribution 
for integrable rectangular resonators was found to be distorted near the ori- 
gin. This effect was later explained []1T| as consequence of the perturbation 



due to the measuring antenna. 

A complementary — and more general — point of view is to regard the 
resonator with attached antennas as a scattering system and to study its 
transport properties; this approach is also natural in view of the mentioned 
analogy with conductance of quantum dots to which quantum wire leads 
are attached. Solving the scattering problem in the described geometry is 
not easy and one naturally looks for possible model simplifications. Since 
the antennas are thin one can try to describe them as lines coupled to the 
planar resonator in which the field obeys the Laplace equation with Dirichlet 
boundary conditions. 

At a glance the different dimensionality of the configuration space parts 
adds complications to solution of the corresponding wave equation. How- 



ever, a standard technique based on self-adjoint extensions [|12l allows us to 
reduce the task to matching solutions in different parts. We shall describe it 
below restricting ourselves to simplest case with one antenna; a more general 



situation will be discussed in a forthcoming paper |T3|. 

Let us begin with a halfiine attached to a plane; we want to couple the 
corresponding free Hamiltonians in a self-adjoint way. In physical terms this 
requirement means the current conservation. It will be ensured by boundary 
conditions (^ given below which are local, and therefore also applicable when 
the halfiine is replaced by a line segment and/or the plane by any nonempty 
planar region. The halfiine-plus-plane system corresponds to the Hilbert 
space L^(iR_) © L^(iR^) , i.e., the field is described by pairs := (^M of 
square integrable functions on the two parts of the configuration space. For 
the sake of simplicity, and with the quantum mechanical analogy in mind, 
we shall speak about wavefunctions. 

The system has a rotational symmetry, so the partial-wave decompo- 



sition may be employed. In particular, the s-wave component $2(^, V^) = 
(27r)^^/^02(^) is independent of the azimuthal angle. The Hamiltonian acts 
on the wavefunction components as the Laplacian, 



H4> = 7' . (1) 




The self-adjointness means a proper matching of the components at the 
junction. A detailed discussion of this problem involving a full classification of 
admissible boundary conditions was presented in |T^. They can be expressed 
in terms of generalized boundary values for wavefunctions in the plane fl^ 

Lo(<l>) := lim-^^, 
r^o inr 

Li(<l>) := lim [$(f)-Lo($)lnr] . (2) 

Typical boundary conditions then read 

0;(o-) = A0i(o-) + 5Lo(<i'2), 

Li($2) = C0i(O-)+DLo($2); (3) 

there are additional lower-dimensional families corresponding to the cases 



where (y) becomes singular |T^. Let us remark that using the appropriate 
regularized boundary values, one can couple in a similar way a halfline to 
IB?, and more generally, two regions whose dimensions do not differ by more 
than two. 

The coefficients A, _B, C, D were in the last named paper expressed in 
terms of elements of the unitary matrix relating the deficiency subspaces. In 
this way they depend on four real parameters; the number reduces to three 
if we demand the Hamiltonian to be time-reversal invariant. A drawback of 
these parametrizations is that they do not show explicitly the ranges for the 
values of the coefficients. 

An alternative approach has been advocated recently in ||I5[. It is based 
on expressing the wavefunction in the plane as $2 = ttGo(-, 0; \/\q) + m(-) , 
where Go is the free Green's function for some energy value Aq which does 
not belong to the spectrum of H and m is a function regular at the con- 
nection point X = . The boundary conditions then tie a and 'u(O) with 



0i(O— ), 0']^(O— ) . In general, one can obtain in this way the same set of ad- 
missible Hamiltonians as above, however, in the mentioned paper a particular 
subset corresponding to A = was discussed and used. 

A more straightforward way is to compute the boundary form to Hq , the 
adjoint to Hq obtained by restriction of a "decoupled" operator, say, the one 
with i? = C = , to functions which vanish at the vicinity of the junction 
[0. Since the action of Hq is given by the same differential expression (|^) 
we find after a simple integration by parts 

(0,ffo»-(^oV,^) = 0'i(O)V^i(O)-0i(O)^;(O) 

+ hm e (He)i^'i{e) - ^',{e)Me)) • (4) 

We know from |[T^ that only the s-wave component in the plane can be 
coupled nontrivially, hence we may suppose that $2(2^) is independent of 
the azimuthal angle and 

02(e) = v^ [^o('^2) Ine + i:i($2) + 0{6) ] . 

Using this asymptotic behavior we can rewrite the the last term on the rhs 
of d) as 

it is then clear that the boundary form is zero under the conditions (P) with 

A, D eM and B = 2tiC . (5) 

If the Hamiltonian should be time-reversal invariant, the boundary condi- 
tions must not change when passing to complex conjugated wavef unctions; 
this happens if the off-diagonal coefficients B, C are real. 

To proceed further we have specify the coefficients in @. The choice 
needs a physical motivation. Since we are interested in transport through 
such junctions, we shall use a comparison of low-energy scattering proper- 
ties with those of a more realistic antenna. The reflection and transmission 
amplitudes for a plane wave e^^^ approaching the junction on the halfline 
were computed in |T2[: we have 



Kfc) = -^, i(fc) = ^-, (6) 



where 



P± := {A±ik) 



2i ( , fc' 

1 + - .-DH-ln- 



TT 



and 7 = 0.5772 ... is the Euler's constant. 

Consider now a semiinfinite cyhndrical tube of radius a attached to the 
plane; we suppose the latter has the central circular area of radius a removed 
so that the whole surface is topologically equivalent to an infinite cylinder. 
The scattering on such a surface is certainly affected by the sharp edge at 
the interface of the two parts but we disregard this fact. 

We have a rotational symmetry again, so each partial wave may be treated 
separately. Suppose that the "longitudinal" component of the incident wave- 
function is a plane wave of momentum k . For an orbital quantum number 
£ one has to match smoothly the corresponding solutions 



fix) :-- 



^ikx _|_ ^^g ikx _ _ _ X < 



^taHf\kr) ... r>a 



We abuse here notation on the Ihs and employ x as a common denomination 
of the longitudinal variable on the cylinder and the radial variable in the 
plane. The matching conditions are solved easily giving 



with 



ra{k) = -'^^ Uk)=uf-^{viy' (7) 



VI := (1 ± 2ika)Hf\ka) + 2ka (HP)\ka) ; 

the Wronskian relation W{Jy{z),Yi,{z)) = 2/7iz implies that the unitarity 
requirement, |ra(A;)p+ |ta(fc)P = 1 , is satisfied. 

It follows from the asymptotic properties of Bessel functions |TB[ that 
|^a(fc)P ~ {kaY^~^ holds for £ > 1 , so the transmission probability vanishes 
fast as fc ^ for higher partial waves. On the other hand. 



Substituting into (|^ we find that these aniphtudes have for fca ^ 1 
same leading behavior as (BI) provided we put 



the 



A:= — 

2a 



D:-- 



Ina . 



BC := - 



(8) 



owing to (I 
satisfied if 



the last relation is in the case of "real" boundary conditions 



(9) 



B = 2-kC = \ — 




In this way we have found values of the coefficients for which the coupling 
@ is able to model the scattering behavior of a real antenna as long as the 
wavelengths involved are much larger than the antenna radius. A similar 
effect has been noticed recently for point interactions in dimension two or 
three [^ : long enough waves feel only the size of an obstacle or a drain. It 
is also worth noting that the coefficient A is always nonzero by (P), so the 
coupling used in Ref.|T^ does not belong to this class. 

Now we are ready to describe our model. We consider a resonator in 
the form of a planar region M to which an antenna is attached at a point 
xo G M ; we suppose that at the junctions the wavefunctions are coupled 
by boundary conditions (H) with the parameter values (||), (P), and that the 
resonator has hard walls which we model by the Dirichlet condition, u{x) = 
at the boundary of M . 

We have to match the plane wave combination e*^^' + re~ 
antenna halfline iR_ with the internal solution u{x) := bG{x, Xi 
G{-, ■; k) is the Green's function of the Dirichlet Laplacian; we assume, of 
course, that k'^ equals to no eigenvalue of this operator. To make use of (|^) 
we have to know the generalized boundary values Li := Li{u; xq) ; they are 



; /c) , where 



— , Li = b^{xo; k) 

ATI 



where 



^{xQ]k):= Jim 

a;— >xo 



In la; — xqI 



27r 



G(x, f o; k) + 

Matching then the solutions we find 

TiZ{k){l-2ika) -1 



r{k) 



'KZ{k){l + 2ika)-l 



(10) 
(11) 

(12) 
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with 



Zik) := axo;k) 



In a 

~2^ 



(13) 



In order to make use of these formulae we have to know the last named 
quantity. If the region M is compact, the corresponding Dirichlet Laplacian 
has a purely discrete spectrum. We denote by A„ amd 0„ , respectively, the 
corresponding eigenvalues and eigenfunctions; without loss of generality the 
latter can be chosen real. This allows us to express the Green's function. 



G{xi,X2]k) = Y^ 



(pn{Xi)(j)n{x2) 



(14) 



n=l 



It diverges, of course, as xi -^ X2] we need to know the regularized value 
(pAl) . A semiclassical argument [|l^. Sec. XIII. 15] makes it possible to asses 



the rate of the divergence: we have A^ ~ 47rn|M| and 



(f)n ^ \M\- 



as ra ^ oo , where \M\ is the area of the resonator. This inspires us to 
employ the identity 



G{xq + i/efi, xq] k) + 

oo ( 



27r 
fo + y/en)(t)n{xo) 



ra=l 



/Krt 1^ 



Aim 



where n is an arbitrary unit vector. The two series diverge as e 
their difference remains finite: we have 



^{xQ]k) = J2 



n=\ 



|0n(^o)P 
An-fc2 



47m 



but 



(15) 



Let us stress that the exchange of the limit and summation should be taken 
with a caution, since the series is not uniformly convergent. Replacing, e.g.^ 
y/e by Cy/e we get an additive factor. However, the correct ^{xq; k) should 
be independent of the boundary at large negative energies. 



lim 



^(fo; ii^) ~ ^y^2^^ 



0, 



by Ref.[|l^]. For a specific M as the one considered below, it is straightfor- 
ward to check that this is true for (p!5|). 
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We are especially interested in the situation where M is a rectangle 
[0, ci] X [0, C2] , for which we have 

(i)nm{x,y) = ^ sin(n— x) sm{m—y) , (16) 

^/c^C2 Ci C2 

the formulae ([T2|) and (p!5D-(p!7|) yield a complete solution of the scattering 
problem in this setting. 

In particular, we want to find the resonance spacing distribution which 
can be compared with the result found experimentally in [^; we refer to 
this paper for a more detailed description of the experimental arrangement. 
By ([T^ ) the resonances are given by complex zeros of the denominator, i.e., 
by solutions of the algebraic equation 

e(fo,A;) = ^+ I (18) 

where a is the radius of the antenna opening. We have solved this equation 
numerically and evaluated the resonance spacing distribution. In order to 
mimick the experimental setting we have evaluated the resonance spaing 
distribution for several rectangular resonators with Ci and C2 ranging from 
20 to 50 cm. The coupling point Xq has been chosen randomly for each 
billiard and the diameter of the antena was equal to 1 mm. The obtained 
data has been combined in order to enlarge the statistics. Since the present 
model was derived under the assumption fca -C 1 we used only the data 
corresponding to frequencies below 10 GHz for which ka < 0.1 . In addition, 
it is known that about 7% of resonances are overlooked in the experiment; 
to take the missing levels into account we have deleted the same amount of 
randomly chosen resonances from the family under consideration. 

The result is plotted on the Figure 1. The agreement between the model 
calculations and the experiment data is convincing. We recall that our argu- 
ment involves no free parameters because the "coupling strength" between 
resonator and the antenna is fully determined by the size of the latter. On the 
other hand, the long-wave condition relative to the antenna radius, ka <^ 1 , 
is important: the agreement worsens if resonances above 10 GHz are taken 
into account. 



In conclusion, we have described a method to treat the scattering prob- 
lem in systems with configuration space composed of parts of a different 
dimensionality. It has been applied to a flat microwave resonator with a thin 
antenna; the calculated resonance spacing distribution agrees in the long- 
wave regime with the experimental data. 
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Figure captions 

Figure 1. The numerically evaluated resonance spacing distribution for 
the rectangular resonator (full line) in comparison with the experimen- 
tal data taken from |ll| (bins). In the inset we sketch the geometric 



arrangement of the experiment. 
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